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ABSTRACT 

We present a physical model for origin of the cosmic diffuse infrared background 
(CDIRB). By utilizing the observed stellar mass function and its evolution as input to a semi- 
empirical model of galaxy formation, we isolate the physics driving diffuse IR emission. The 
model includes contributions from three primary sources of IR emission: steady-state star for- 
mation owing to isolated disk galaxies, interaction-driven bursts of star formation owing to 
close encounters and mergers, and obscured active galactic nuclei (AGN). We find that most 
of the CDIRB is produced by equal contributions from objects at z ~ 0.5 — 1 and z ^ 1, as 
suggested by recent observations. Of those sources, the vast majority of the emission origi- 
nates in systems with low to moderate IR luminosities (Lm <J 10 12 L Q ); the most luminous 
objects contribute significant flux only at high-redshifts (z ;> 2). All star formation in ongo- 
ing mergers accounts for <J 10% of the total at all wavelengths and redshifts, while emission 
directly attributable to the interaction-driven burst itself accounts for < 5%. We furthermore 
find that obscured AGN contribute <^ 1 — 2% of the CDIRB at all wavelengths and redshifts, 
with a strong upper limit of less than 4% of the total emission. Finally, since electron-positron 
pair production interactions with the CDIRB represent the primary source of opacity to very 
high energy (VHE: E 1 <; 1 TeV) 7-rays, the model provides predictions for the optical depth 
of the Universe to the most energetic photons. We find that these predictions agree with ob- 
servations of high-energy cutoffs at ~ TeV energies in nearby blazars, and suggest that while 
the Universe is extremely optically thick at J; 10 TeV, the next generation of VHE 7-ray 
telescopes can reasonably expect detections from out to ~ 50 — 150 Mpc. 

Key words: galaxies: starburst - galaxies : active - galaxies: formation - cosmology: diffuse 
radiation - infrared: galaxies - gamma-rays: theory 




therein). However, in regions with significant dust opacity radi- 
ation generated by stars and AGN i s reprocessed by dust into 
the IR as thermal emission (see, e.g. Soife r & Neugebauedfl99ll : 



Diffuse extragalactic background light (EBL) represents the sum 
total of all the photons produced by luminous matter over the 
lifetime of the Universe. In the ultraviolet, optical, and near- 
infrared (IR) the EBL is directly attributable to star for mation 
and active galactic nucleus (AGN) activity (for rev iews, see lTvsoJ 



1 INTRODUCTION 



Sanders & Mirabel|[l996h . leading to the expectation of a signifi- 
cant EBL component at A ~ 10 — 1000/im. Measurements of this 
Cosmic Diffuse IR Backround (CDIRB) were notoriously difficult 
to obtain, owing primarily to the presence of significant zodiacal 
and galactic foregrounds and the lack of access to this region of 
the spectrum from the ground. The launch of the Cosmic Back - 
ground Explorer (COBE: for an overview, see lBoggess et al.ll992h 
finally revealed a CDIRB comparable in brightness to the optical 
EBL, providing a complete censu s of obscured star formation and 
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Because the emission mechanisms that generate the CDIRB 
are intimately connected to galaxy formation and evolution, 
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it provides a powerful observational constraint on models 
jPartridge & Peebles! [l967). A particularly popular technique 
for modeling the CDIRB ha s been backwards evolution (e.g. , 
Rowan-Robinspnl 1200 ll. 120091; lLagache et allliool |2004 12005 ' 



XuetalJl200~ ; iFranceschini et alj|2008l ; iFinke & Razzaquell2009 
in which a parameterized fit to the evolution of the IR lu- 
minosity function - using low-redshift observations as a base- 
line - to reproduce a number of observables, including number 
counts and the CDIRB itself. A compl ementary approach has 
been to use semi- analytic mo dels (SAMs: ICole et alj|l994 [2000: 
ISomerville & PrimaclJ [l999T) . in which the hierarchical growth 
of structure is tracked by N-body cosmological simulations and 
baryonic physics - e.g., star formation, radiative cooling, gas ac- 
cretion, AGN activity, etc. - are implemented as simple ana- 
lytic prescriptions, tuned to match a predetermined set of (usu- 
ally local) observational constraints. Mock galaxy catalogs gen- 
erated by SAMs can then be combined with template spectral 
energy distributions (SEDs) to construct predictions for number 
counts at a given waveleng t h and t he CDIRB more generally (e.g 
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Primack et all Il999l 12005. 20081: iDevriendt & Guiderdonil |200C|: 



Baug h et al .l2005l : ISwinbank et alteOolT Finally, c osmic chemica l 



evolution (CCE) techniques (e.g.. iPei & Fal l 1995; Pei et al. 1999) 
model the EBL from the optical through IR by solving for the self- 
consistent evolution of globally averaged quantities - such as gas 
depletion and star formation - i n an analogou s manner to galactic 
chemical evolution models (e.g.. lTinslevlll98rj|) . 

Each of these approaches has clear advantages and disadvan- 
tages. Backwards evolution models, while robust and successful 
in reproducing the observations, must assume simple evolution- 
ary scalings. Therefore, they contain no physics, and cannot sep- 
arate out the relative contributions to the CDIRB from different 
emission mechanisms. SAMs, by contrast, can determine the rel- 
ative importance of, e.g. AGN versus star formation, but rely on 
simple analytic implementations that introduce parameter degen- 
eracies and may not include all of the complexities in the rele- 
vant physics - though some have sought to address this shortcom- 
ing by incorporating the results of hydrodynamical simulations of 
merge r driven starbursts an d AGN activity into the SAM frame- 
work dSomerville et alJ2008h . Furthermore, SAMs have also faced 
difficulty simultaneously matching the o bserved galaxy mass fu nc- 
tion and cosmic star formation history dSomerville et alj[2008l) ; it 
is no generally possible to separate these well-known issues from 
other model elements in driving their predictions for the CDIRB. 
Finally, while CCEs do not require a detailed implementation of 
these complicated baryonic physics, their global nature means they 
have little to say about specific populations of objects. 

An alternative approach, outlined in detail by iHopkins et al.l 
(2008b, a) attempts a semi-empirical model of galaxy formation in 
the context of a merger-driven cosmic cycle: while isolated galax- 
ies domina te the overall cosmic energy budget (however, see dis- 
cussion in Hopk ins et al. I l2010al) , major mergers of gas-rich disks 
drive an evolutionary sequence through a period of intense star 
for mation and AGN act i vity, to produce passi ve elliptical galax- 
ies dSanders et ail 1 19881 ; IHopkins et alj|2006bh . The model con- 
nects measurements of the stellar mass function to the dark mat- 
ter halo populat ion through a halo occupation distribution (HOD) 
approach (e.g., Peacock & Smith 2000; Scoccimarro et al. 2001; 
Berlind & Weinberg 2002; Berlind et al. 2003; Kravtso v et al 



2004 ; jZehavi et aiTl2004 IZheng et all I2005L 120091 : iBrown et all 
2008MConrov & Wechslej2 009). and implements baryonic physics 



served quasar luminosity f unction (IHopkins et al.l [2008b), the 
growth of the red sequence dHopkins et al. 200 8^), the evolutio n 
of massive elliptical galaxies with redshift ( IHopkins et al1 l2009d). 
bulge-to- disk rations {H opkins et al. 2009f), and the IR luminosity 
function [Hop kins et alj2010bl) . Semi-empirical modeling also rep- 
resents a fundamental advance over pervious techniques for mod- 
eling the CDIRB; it is predictive, incorporates the relevant physics, 
includes the results of high-resolution hydrodynamical simulations 
where appropriate, while at the same time is constructed to match 
the galaxy mass function at all redshifts. Therefore, we can isolate 
the physics driving the production of the CDIRB without signifi- 
cant degeneracies model components and parameters. 

In this work, we use this semi-empirical framework to con- 
struct a physical model for the origin of the CDIRB. In addition 
to providing an important verification of this particular model of 
galaxy formation more generally, this approach will allow us to 
unfold the contributions from objects as a function of luminosity 
and redshift, as well as by emission mechanism - including star 
formation in isolated disks, merger-driven starbursts, and AGN ac- 
tivity. Furthermore, because electron-positron pair production in- 
teractions (7 + 7 — i> e + + e~) with CDIRB photons is the pri- 
mary source of attenuation of extragalactic, ~ TeV photons (of- 



ten referred to as Very High Energy, or VHE 7-rays;|J elley 



Fazio et al. 1966; Fazio & Stecker 1970; Gould & Schreder 



1966 



1967 



Stecker 1969; Stecker et al. 1992), model predictions for the evo- 
lution of this background with redshift yields a prediction for the 
opacity of the Universe to the most energetic photons. 

This work is organized as follows: in § [2] we outline our 
methodology, in § [3] we summarize the primary sources of un- 
certainty in the model predictions, in § [4] we verify that the 
model is consistent will the relevant observational constraints, in 
§ [5] we investigate the result of changing some of the model as- 
sumptions, in § [6] we present predictions for the relative impor- 
tant of different populations of objects in generating the CDIRB, 
in § |7] we make predictions for the 7-ray opacity of the Uni- 
verse, and in § [8] we conclude. Throughout this work, we will 
make frequent reference to the total IR luminosity which is de- 
fined according t o convention as the integ rated luminosity from 
8-1000 /im (e.g.. ISanders & Mira bel 1996). Furth ermore, we as- 
sume t he most recent cosmological parameters from lKomatsu et all 
d2010l) : (Q m ,tt A ,h) = (0.26,0.74,0.71), the WMAP 7-year 
mean. However, it is important to note that varying these within 
a reasonable range has no noticeable effect to our results. 



by incorporating the results of high-resolution hydrodynamical 
simulations. This approach has successfully reproduced the ob- 



2 METHODOLOGY 

In this work, we make extensi ve use of the semi-empi rical model of 
galaxy ev olution presented by IHopkins et al . ( 2008b a), and subse- 
quently bv lHopkins et al. ( 2009d, 2010b, 2009a 1); i n partic ular, we 
utilize the model IRLFs presented bv lHopkins et alJd2010bl) . These 
predictions, which are consistent with the observed galaxy mass 
function, IRLF, and cosmic star formation history from z ~ — 4, 
include contributions from the three dominant sources of emission: 
steady-state star formation, merger-driven starbursts, and obscured 
AGN activity. While we refer the reader to the above-mentioned 
references for details, here were briefly summarize the principal 
model components and methods and how they are used to construct 
predictions for the CDIRB. 
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2.1 Steady-state star formation 



Given the obser ved 
Bell et all 120031: 



at redshift z (e.g. . 



Perez-Gonzalez et al 



stellar mass funct ion 

20061. lArnouts et alj 120071 : 
we populate 



Fontana et al. _ 

20081 : lllbert et alj l20ld) 



dark matte r halos with galax i es fo llowing the methodology 
outlined by IConrov & Wechslerl J2009h . These galaxies are then 
assigned properties consistent with observations of systems at that 
redshift with comparable stellar masses. 

The gas mass is estimated from observations of the gas 
mass fractions (/ 9 ) of star-for ming disks, and is a strong func- 
tion of mass and redshi f t (e.g., iBell & de Jong 2001; Kannappan 
20041: iMcGaughl 120051 ; IShaplev et alj l2005UDaddi et alj 1200 J 



Tacconi ct al. 



2010) 



It has been shown ( Hopkins et alj 12008b 



2009d li ) that this evolution is well-described by a simple functional 



form, motivated by cosmological simulations (Kere s et"aT]|2005l 
l2009h : 
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where fw(z) is the fractional look-back time from the present day 
to z (fit — > 1 as 2 — s> oo). Motivated by the observations listed 
above, we furthermore assume a constant scatter of ~ 0.25 dex 
independent of stellar mass. 

The size of a given disk galaxy is assigne d consistent with 
the observed size-mass relation ( Shen et al. 2003) and its evolution 
(e.g.. iTruiillo et alj|2004l ; iFerguson et alj|2004l : iToft et alj|2007l) . 
ISomerville et alj ( 120081) found that this too can be represented by a 
simple analytic fitting form: 

/ il/f \ 0.25 

R e (M ir \z = 0) = Ro ~ 5.28 kpc x ( - 



Mi, 



R^M* 



Ro(l 



(2) 



Again, motivated by observations we assume an intrinsic scatter of 
0.2 dex. 

Finally, given the gas mass and size of a disk galaxy, the 
rate of steady-state star formation - i.e., quiescent, or owing pri- 
ma rily to internal gas dy namics rather than external factors (see 
also lHopkins et alj|2010al) - can be estimated u sing the Kennicutt- 
Schmidt Law ( Schmid il959l ; lKennicuttll 1998bl) : 



S s/ = (1.3 x 10" 



©yr 



1 kpc" 



(3) 



where E fl is the disk-averaged gas surface density in units of 
Mq kpc -2 , £ s / is the disk-averaged star formation rate surface 
density, r?x = 1.4 fo r consistency with observations (see also 
lKenmcutd[l989l,[ 1998a), and the normal zation has been corrected 
2003T) stellar initial mass function (IMF). 



Chabrier 



assuming a 

Given the assumed IMF, we can convert this star formation rate 
to the total IR l uminosity using the known linear scaling (e.g., 
lKennicuttlll998al) . 



2.2 Interaction-driven starbursts 

In addition to steady-state star formation, close encounters between 
and mergers of gas-rich galaxies can concentrate large quantities of 
gas in the nucleus, l eading to a dramatic enhancement in the total 
star formation rate iHernouist & Katzl 1 19891 : iBarnes & Hernquisl 
Il99ll : iMihos & Hernquistlll994l Il99fjf) . This owes principally to 
the effects of the time- varying gravitational potential during close 
passages: an intrinsic lag in the induced stellar relative to the 




Rest Frame Wavelength [/xm] 



Figure 1. Examples of the template SEDs use to connect the model IR 
luminosity functions to observable quantities. We show the templates for 
star-forming systems (black) and AGN (red) at four different luminosities: 
log(L /fi /L Q ) = 10 (dotted line), 11 (dashed line; a LIRG), 12 (solid 
line; a ULIRG), and 13 (dot-dashed line; a HyLIRG). The temp lates for 
star-forming s ystems are from the library o f lCharv & Elbaj 120011) . and the 
AGN are from Sieben morgen et all l2004bl and references therein) 



gaseous bars drains the gas of its angular momentum within a 
charac teristic radius, sending it on a free-fall trajector y to the nu- 
cleus (Barnes & Hernquist 1996; Hopkins et al. 2009c). The result- 
ing nuclear gas concentration fuels a massive starburst that can 
enhance the global starburst by more than an order of magnitude 
dPi Matteo etafll200l 120081 ; ICox et alj|2008h , and yields charac- 
teristic central mas s concentrations that have bee n observed in local 
elliptical galaxies ( Hopk ins et all2008cl l2009el lbh . 

Mergers in our model are tracked via halo-halo merger 
rates estimated from cosmological N-body simulations (e.g., 
IStewart et alj2008l ; lFakhouri & Mall200cll2009l) . We then estimate 
the galaxy-galaxy merger rate at z as a function of stellar mass, 
redshift, and mass ratio {p = Mi /M2 < 1) by first correcti ng for 
the dynamical friction time delay (Bovlan-Kolch fnet al.l2008l) . and 
then convolving with the galaxy population using the HOD model. 
Though the exact star formation history of an individual encounter 
is impossible to determine from our model, the IRLF requires only 
a knowledge of the distribution of times spent at a given star forma- 
tion ra te, which can be estimated using an ensembl e of simulations 
(e.g jHopkins et alj|20063) . lHopkins et alj ( feOlObl) found that this 
can be represented by the fitting function 



dt 



dloeM* 



t a b exp - 



M ab /t a 



(4) 



where M B f, is the total mass participating in the starburst , and 
t a b ~ 100 Myr is the characteristic timescale dCox et al.l2 008). It is 
important to note that tracking the gas content of galaxies is a criti- 
cal component of the model. Owing to simple gravitational physics, 
the burst properties can be represented in simple analytic form (in 
agreement with the results of hydrodynamical simulations); the to- 
tal gas mass that is stripped of its angular momentum and par- 
ticipates in the nuclea r starburst scales as M a b ~ (1 — f g )fj,M g 
( Hop kins et ai] |2009c ) . In the context of the model, star formation 
owing to the interaction-driven burst itself is independent of the gas 
that retains its angular momentum and continues steady-state star 
formation in the disk - a distinction that will become important 
when characterizing the relative importance of these two modes of 
star formation (e.g., in § 16.21 1. 
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2.3 Obscured AGN 

In addition to a nuclear starburst, it has been shown that en- 
counters between gas-rich galaxies will excite AGN activity 
dDiMatteoetalj2 005; Hopk ins et al.l2005bl l5). Furthermore, more 
recent analyses of hydrodynamical simulations have shown that 
the peak luminosity of these AGN, corresponding to the Edding- 
ton limit of the final black hole mass, is proportional to the 
binding energy of t he bulge in the remnant dHopkins"et ai] |2007l ; 
lYounger et al.l2008h . This bulge mass consists of both the gas mass 
participating in the burst (M s &) and the fraction of the stellar mass 
which undergoes violent relaxation (Lvnden-Bell 1967). This frac- 
tion is proportional to fiM* , again owing to scalings derived from 
simple gravitational dynamic s and in agreement wi t h the results of 
hydrody namical simulations (Hopkins et al. 2009c). Hopki n~et al.l 
(2010b) find that the relationship between the peak AGN luminos- 
ity and these quantities can be conveniently approximated by: 



.AGN 



4.6 x 10 11 L Q x (1 + zf 



M sb + jj,M± 
10 10 M© 



(5) 



where the redshift evolution reflects the effects of higher 
as fractions and s maller scale-lengths in high-redshift disks 
Hopkins et ai]|2007l) . The distribution of time spent above a given 
luminosity has been studied extensively in a series of papers (see 
Hopki ns et al. 2006bl. l2008bl and references therein), and is well- 
described by: 
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Thus, given the merger-rate as a function of redshift, stellar mass, 
and mass ratio we can estimate the peak luminosity and convolve 
with the distribution of time above Ltd to yield an estimate from 
th e bolometric luminos ity function (the whole exercise is presented 
in iHopkins et alj[2010bl) . However, while it is trivial to convert 
a star formation rate to an IR luminosity, the situation for ob- 
scured AGN is rather more complicated. Therefore, to transform 
the bolometric AGN luminosity function to t he IRLF, we ado pt the 
luminosity-dependent luminosity fraction of lGilli et al. (2007) and 
assume this obscured luminosity is reprocessed into the IR. Though 
this is a simplification, we find that it yields simil ar results to either 
using the observed column density distribution dUeda et al.ll2003h 
or the mean X-ray to IR bolometric correction for obscured AGN 
dElvis et al.lll994l ; IZakamska et alj|2004 IPolletta et alfepOd) . 



2.4 Predicting the diffuse background emission 

Given a theoretical IRLF (dn/dlogLiR) and an assumed cosmol- 
ogy, the predicted CDIRB is: 



/„ (CDIRB) 



dz I d\ogLiR 



S u (Lir,z) x 

dn dV 
dXogLm dzdVL 



(8) 



where S v (Lir,z) is the flux density at frequency v in the ob- 
served frame for an object with total IR luminosity Lir at redshift 
z, and dV/dzdQ, is the comoving volume element. For this work, 
we derive S(Lir,z) for star formation dominated systems (both 
steady-state and merger-driven star formation) by assigning them 



templates from the library of ICharv & Elbaj d200ll), whi ch them- 
selves include models first presented in lDale et al. [ feOOlh models. 
This template library has been found to reproduce many of the scal- 
ing relations observed in local IR-luminous galaxies. For clarity, 
examples of these templates for three different luminosities are pro- 
vided in Figure [Tj 

The choice of SED templates for obscured AGN represents 
a somewhat more complicated exercise. While a simple empirical 
template - e.g., the local ULIRG Mrk 231, which owes much of its 
luminosity to a deeply obscured AGN - can be used, this approach 
does not account for potential dependence of the effective dust tem- 
perature on t he luminosity of the system, as is the case for the star- 
burst models dCharv & Elbad200ll ; lDale et al.l200ll : lDale & Helot] 
l2002l : lLagac"he' et alj|20 03). Therefore, we choose t o utilize a well- 
studied suite of models for heavily obscured AGN Jsiebenmorgenl 
199 lL Il993l; ISiebenmorgen & Kruegell 1 19921 ; ISiebenmorgen et all 
1 992II2OO1I , 2004a b Fl which are parameterized by the outer radius 
of the dust clouds (R), the total emergent IR luminosity (Lir), and 
the visual extinction (Av) - or alternati vely the hydrogen colum n 
Nh/Av » 1.8 x 10 21 cm~ 2 mag" 1 dPredehl & Schmittlll995h . 
In particular. ISiebenmorgen etaL (2004b) find that optical Type II 
AGN are best fit by A v « 16 - 32 mag -> \og(N H / cm~ 2 ) » 
22.5 — 22.8, and R ~ 1 — 2 kpc. These columns are also consistent 
with the observed distributi on of column den sities in AGN, inferred 
from X-ray spectroscopy dUeda et al]|2003h . Given these observa- 
tional constraints, for most of this work we assume R = 1 kpc and 
Ay = 16, or log(A r if/cm -2 ) ~ 22.5; we include examples as a 
function of luminosity in FigureQ] 



3 SOURCES OF UNCERTAINTY 

There are three primary sources of uncertaint y in the model results : 
(1) the stellar mass function (also shown bv lHopkins et alj EzOlOb. 
this dominates the uncertainty in the predicted IRFLs), (2) the tem- 
plate SEDs of star formation dominated objects, and (3) the tem- 
plate SEDs of obscured AGN. The first owes to the dynamic range 
of observational estimates for the stellar mass function at different 



redsh i fts (e. g., Bell et al. 2003; Bundv et alj 120051; iFo ntana et al 
20061 ; Ipjrez-Gonzalez et alj|2008l ; I Arnouts et al.ll2007l; lllbert et al 
20 id) w hich is the primary model input. As in [Hopkins et al 
d2010bh . we represent this uncertainty by characterizing the 
full dynamic range spanned by the observations (see § 2.1 of 
Hopkins et al. 20 1 Obi for a detailed description), which essentially 
defines the upper and lower bounds for the model predictions. This 
uncertainty naturally increases with luminosity and redshift owing 
to the lack of large, robust samples of objects which yield a larger 
effective dispersion among different observational techniques. In 
the context of the CDIRB, this yields a ~0. 15 dex uncertainty 
in intensity (/„) at X b s ~ 100 /im , increasing to ~0.2 dex at 
longer wavelengths; this owes to larger systematic uncertainties in 
the IRLF at the high luminosities and redshifts that produce this 
part of the CDIRB (see § I6.lt . It is also important to note that an 
a pparent spike in integ ral quantities near z ~ 2 (e.g., Figs. 4 & 8 
in Hop kins etalHo iOb) is an artifact of observational estimates for 
the stellar mass function at these redshifts and not a physical effect 
owing to the model itself. 

While we find that the observed mass function represents the 
primary source of uncertainty in model results, for this work we 



1 Available at htt p : / /www . eso ■ org/~rsiehenm/agn_models/| 
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Figure 2. The predicted CDIRB from the model (solid red line), includin g the uncerta i nty in normalization (red shaded region) as compared to observations: 
the black shaded region shows the fit and normalization uncertainty from Fixsen et al. 1 1998) to data obtained by the FIRAS experiment on board COBE, the 
filled diamonds upper limits from 7-ray absorption (see §[7]for more details : lMazin & Raudl2007n . and and the remaining observations are listed in Tablel2.3l 



Table 1. Observational Constraints on the CDIRB 



Wavelength 


Symbol 


Instrument" 


Reference 


1-100 /im 


Filled diamond 


IACTs 


Mazin & Raue (2007) 


15 (im 


Lower limit 


ISOPHOT 


Chary &Elbaz (2001) 


24 fim 


Open diamond 


MIPS 


Papovich et al. (2004) 


60 fim 


Upper limit 


IRAS 


Miville-Deschenes et al. (2002) 


60, 100 fj,m 


Inverted triangle 


DIRBE 


Finkbeiner et al. (2000) 


65, 90, 140, 160 fj,m 


Pentagon 


AKARI 


Matsuura et al. (20 1 0) 


70, 160 /im 


Asterisk 


MIPS 


Dole et al. (2006) 


100-1000 fj,m 


Black region 


FIRAS 


Fixsen et al. (1998) 


140, 240 fim 


Square 


DIRBE 


Hauseretal. (1998) 


150-180 lira 


Cross 


ISOPHOT 


Juvela et al. (2009) 


170 /xm 


Star 


ISOPHOT 


Laeache & Puget (2000) 


250, 350, 500 fim 


Circle 


BLAST 


Patanchon et al. (2009) 


850 fim 


Filled Triangle 


SCUBA 


Knudsen et al. (2008) 



a Instruments include: Upper limits from imaging atmostpheric Cherenkov telescopes (IACTs: Finlev & The VERITAS Collaboration 2001; Daum et al. 
ll997tlHintorJl2004l ; ICortina et alj2005l) . the Multiband Ima ging Photometer for SIR TF (MIPS: lRieke et alj2004l) on board the Spitzer Space Telescope 
(Werner et al. 2004), the Infrared Astronomy Satellite (IRAS- lNeueebauer et aljl984h . the Far-Infrared Absolute Spectrophoto meter (FIRAS) and Diffuse 
Infrared Background Ex periment (DIRBE) o n board the Cosmic Background Explorer Satellite (CO BE: for an o verview. seelBoggess et all 19921) . the 
Far-Infrared Surveyor (FIS: lKawada et alj2007h on board the AKARI Sattellite ^Murakami et al.l2007h. ISOPHOT jLemke et aljl996h on board the Infrared 
Space Observatory (ISO), the Balloon-borne Large-Aperture Submillimetre Tele scope (BLAST: IPascale et al.l2008l) . and the Submillimetre Common User 

Bolometer Array (SCUBA: [Holland et all 19991) . 



also must consider our lack of knowledge of the template SEDs. 
This can be done in a similar manner to the mass function uncer- 
tainty, by considering alternative suites of template SEDs and/or 
considering limiting cases. For star formation dominated objects 
- including both the steady-state and merger-drive n modes - we 
compa re the results of using the template l ibrary of lLagache et al.l 
( 120031) as opposed to our fiducial choice o f lCharv&ElbazlfeOOil) . 
This yields a ~0.05 dex normalization uncertainty independent of 
wavelength. We add this in quadrature to our model predictions 



in Figure [2] but note that the model uncertainty is still dominated 
by the dynamic range in stellar mass function estimates. Obscured 
AGN again represent a more complicated case, and we consider 
two methods for determining the systematic uncertainty associated 
with their contribution to the CDIRB: a limiting case in which the 
full bolometric luminosity is reprocessed into the IR, and the effects 
of varying the templated SED model parameters (see § |6.3t . 
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Figure 3. The cumulative fraction of the CDIRB resolved as a function of flux density. The left pan el shows model predi ctions for 24 fim (blue), 70 /im 
(re d), and 160 fim (b lack) as com pared to observati ons: the black ci rcle is from [ Fray er et alj ^20063) , the tri angles [FraveretalJ l2006bl) , the inverted trian- 
gle [Dojeetal] (20061) . the squares IPoIe et alj 120041) . the filled stars iBethermin et alj j2010|) , and the asterisk IPoie et alj feOOll) . The right panel shows the 
BLAST/Her schel bands at 2 50 (im (black), 350 fim (blue), and 500 (red) bands as compared to observations: the vertical lines indicate the error range 
from lPatanchon et alj <2009l) . Finally, the bottom panel sho ws the 850 fim (SCUBA2) and 1 1 00 (Az TEC) bands as compared to observations: the square, 
circle, and triangle are estimates from the num ber counts of IWeiB et aTl l2009l) , ICowie et al] j2002l) . and Knudsen et al. 1 2008) respectively, and the star is an 
estimate from the number counts presented by Auster mann et alj 12010ft . 



4 COMPARISON TO OBSERVATIONS 

The primary observational constraint is clearly the observed inten- 
sity (/„) as a function of wavelength of the CDIRB itself. In Fig- 
ure [2] we present the model predictions, including the full range 
of uncertainty, as compared to a range of observations (see Ta- 
ble l2.3l for references). We find that the model predictions are con- 
sistent with direct observations from the COBE satellite at w 100 — 
lOOO^im, as well as estimates from number counts in the same 
wavelength ranges from instruments like IRAS, MIPS, BLAST, 
and SCUBA. While many different observations seem to agree in 
this wavelength range, the mid-IR from ~ 10 — 100 fan has been 
a source of considerably more controversy owing to the difficulty 
of effectively subtracting galactic and zodiacal foregrounds and the 
general lack of access to this wavelength range from the ground. 
In particular, at 60 fim upper limits from th e power-spectrum of 
IRAS maps dMiville-Deschenes et ail |2002|) and the TeV opac- 
ity o f intergalactic space owing to electron-positron pair produc- 
tion (ISteckeretalJll992k iDwek & Krennrichll2005l ; iMazin & Raul 
120071 and references therein) a re potentially inconsiste nt with a re- 
analysis of the DIRBE data by Finkbeiner et alj (2000) - the latest 
analysis bv lMatsuura et alj fcOlfJ has fairly large systematic uncer- 
tainties, and thus does not favor either scenario. Our model predic- 
tions appear not to strongly rule out this measurement, but clearly 



indicate a preference for a dimmer CDIRB in the mid-IR than they 
suggest. 

Another important test of the model predictions is the relative 
contributions from individual sources as a function of their flux 
density. In Figure[3]we present this comparison to observations for 
24, 70, 160, 250, 350, 500, 850, and 1100 data. It is clear that 
the model can account for the vast majority of the CDIRB; at 24, 
70, 160, and 850 fim we find good agreement with the fraction of 
emission at those wavelength above a given flux density down to w 
10% of the total emission. At the same time, however, we find that 
the model cannot account for the ~ 10 — 15% of the CDIRB owing 
to the brightest sources at X b s J> 250/^m. This is a well-known 
problem, exemplified by the difficulty semi-an alytic models have 
found in reproducing the bright 850 fim counts (Bau gh et alj2005l ; 
ISwinbank et alj |2008). There are a number of potential solutions 
have bee n proposed in the literature, including cosmological gas 
accretion dDave et alj20 09) and a more sophistica ted calculation of 
the IR SED including full radiative transfer (e.g.. lNaravanan et alj 
l2010l) combined with the well-known effects of cosmic variance in 
these rather small-area survey s (see discussion in e,g,. IWeiB et alj 
l2009l ; lAustermann et al.l2oTol) . However, because the model agrees 
with observations over > 90% of the CDIRB emission, we do not 
believe this discrepancy materially affects our results. Modeling the 
bright end of the luminosity function and the number counts more 
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Figure 4. The predicted CDIRB under three different sets of model assump- 
tions: the fiducial model (red), a steeper Kennicutt-Schmidt Law (rjfc ~ 1 .6 
rather than 1.4; blue), and no significant evolution in the gas content of disk 
galaxies with redshift (black). 

generally is, however, clearly a very interesting problem in its own 
right and a full analysis is currently ongoing (C. C. Hayward et al., 
in preparation). 



5 CHANGING THE MODEL ASSUMPTIONS 

That our model is entirely consistent with available observational 
constraints is very encouraging. However it is important to ex- 
amine how its predictions for the shape and normalization of the 
CDIRB are affected by changing our assumptions. For example, 
strong evolution in the gas content of galaxies and the constancy 
of the Kennicutt-Schmidt Law with redshift - both critical compo- 
nents of the model - may themselves drive the model predictions. 
Relaxing them provides insight into the origin of the CDIRB, as 
well as a test of the robustness of our results. It is also important to 
note that while we could just as well vary the model assumptions 
around the sizes and evolution of disks, the trend of obscured frac- 
tion with AG N luminosity a n d redsh ift, and evolution of the IMF 
with redshift, Ho pkins et alj J2010bh found that these were not a 
significant source of uncertainty in their predictions. 

In Figure [4] we present the model predictions for the SED of 
the CDIRB for two scenarios: one in which the gas content of disks 
does not evolve with redshift, and another in which we assume 
a somewhat steeper Kennicutt-Schmidt Law (the effects of both 
of thes e variations on the IRLFs are presented by Hopk ins~et al.l 
l2010bh . Clearly the most dramatic differences arise from ignoring 
the redshift evolution of gas fractions at fixed stellar mass. This 
systematically suppresses the Universal star formation rate at inter- 
mediate to high redshift relative to the fiducial model by lowering 
the gas mass of those disks. Its effect on the predicted CDIRB is 
two-fold: first the normalization - particularly at J> 100/im - is 
substantially lower, reflecting the systematic decrease in Univer- 
sal star formation rate density, and therefore IR luminosity den- 
sity at z > relative to the fiducial model; second, it shifts the 
peak of emission to shorter wavelength, owing to the preferential 
suppression of high-redshift star formation which is produced by 
emission long ward of the peak in a typ ical starburst SED (X b s ~ 
60 - 100Mm; ISanders & Mirabel 19961) . While there is strong ob- 
servation al evidence that high-redshift disks are preferentially gas- 
rich (e.gjBell & de Ionj|200ll;[Kannappadl2004l;[McGaugrif2005l ; 
IShaplev et alj2005l ; lDaddi et al]2009l ; lTacconi et al.l20ld) . this ex- 



Figure 5. The differential contribution to the predicted CDIRB as a function 
of redshift including all objects (black) and separated into three luminosity 
bins: 'normal' (blue), LIRGs (green), and ULIRGs (red). It is also important 
note that the apparent discontinuity near z « 2 is an artifact of the observa- 
tional estimates for the stellar mass function at those redshifts adopted for 
the model, rather than a physical effect. 

ercise demonstrates our model requires a systematic and substantial 
evolution in their gas content to remain consistent with observa- 
tions. Put another way, we can interpret the shape and normaliza- 
tion of the CDIRB as reflecting this increase with redshift of the 
typical gas fraction of galaxies at constant stellar mass. 

Though local star-forming galaxies are well-fit by a Kennicutt- 
Schmidt index of r\K ~ 1-4, there is some observational evidence 
for a steepe r relation at high-redshift and in more extreme envi- 
ronments (Bouc he et al.l2007l ; lBothwell et al.l2010l) . Therefore, we 
consider a model which is otherwise identical to the fiducial set 
of assumptions, but with t\k ~ 1-6 (renormalizing t he relation to 
avoid o verproducing the local star formation rate). Hopki nTet al.l 
( bOlObh found that the primary effect of this change was somewhat 
more star formation in the most massive systems at high-redshift. 
The result is a set of IRLFs that are largely unchanged considering 
the systematic uncertainties in the model, though with a slightly 
higher contributions form the high-luminosity, high-redshift tail. In 
the context of the CDIRB, the results are again similar with a bit 
more emission at A D (, S <; 300/im - a region of the SED which 
we will see in the proceeding section has the most significant con- 
tribution from high-redshift ULIRGs. But again, considering the 
systematic uncertainties, we find that the model predictions are rel- 
atively insensitive to this change in t\k ■ 



6 INTERPRETING THE COSMIC DIFFUSE INFRARED 
BACKGROUND 

6.1 Sources of emission 

Given that our model predicts a CDIRB that is consistent with the 
observational constraints, we can use it to interpret the origin of 
the emission. In particular, we can separate out the relative con- 
tribution from sources at different redshifts and in three different 
luminosity bins: 'normal' systems {Lir < 10 11 L©), luminous 
infrared galaxies (LIRGS: 10 11 L@ < L IR < 10 12 L Q ), and ul- 
tralumionous infrared galaxies (ULIRGs: Lir > 10 12 L©). In Fig- 
ure[5]we show the differential contribution of each to the integrated 
background emission as a function of redshift. The model shows 
that observed intensity of the CDIRB is dominated by sources at 
z ~ 0.5— 1 and consists primarily of 'normal' systems and LIRGs. 
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Figure 6. Model predictions for the fraction of CDIRB flux density as a 
function of wavelength owing to sources at (top to bottom) redshift of z > 
0.5, 1.0, 1.5, and 2.0 including objects of all luminosities. For comparison 
we show BLAST observations of the fraction of CDIRB lig ht contribute d 
by sources at z > 1 estimated via a stacking analysis (Marsd en et alj200^) . 



dMarsden et afll2009l : I5evlin et alj|2009l) . We also find a greater 
contribution from high-redshift systems at longer wavelengths past 
the peak of the CDIRB (A > 200pm); this again makes sense con- 
sidering the peak in emission for IR-luminous systems that dom- 
inate th e local IR luminosity den sity is ~ 60 — 120 fim (see re- 
view bv lSanders & Mirabelll99r3) . and therefore longer wavelength 
emission in the CDIRB must arise primarily from redshifting this 
peak to longer wavelength. 

We perform a similar exercise in Figure [7] breaking down the 
CDIRB into contributions from the three luminosity bins. We find 
that LIRGS contribute ~ 50% of the observed emission at all wave- 
lengths; 'normal' systems contribute ~ 40% at A (, s ~ 20/im de- 
creasing to ~ 30% at X t s ~ 1000pm while the contribution from 
ULIRGs increases from ~ 5% to ~ 20% over the same interval. 
This owes to the strong evolution in the IRLF of ULRGs, yielding 
a much more important contribution at z > 1 than at low-redshift 
dLe Floc'h et alj|2005l; iMagnelli etaT]|2009h and thus contributing 
a much larger fraction of the CDIRB at wavelengths that are domi- 
nated by emission from high-redshift objects. 
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Figure 7. Model predictions for the fraction of CDIRB flux density as a 
function of wavelength owing to sources at in three different luminosity 
bins: 'normal systems' (blue), LIRGS (green), and ULRIGs (red). 



In general, since dI IR /dz ~ p,(z)(dV/dz)/D 2 L (z) - where p* 
is the cosmic star formation rate density, dV/dz is the comoving 
volume element, and Dl is the luminosity distance - this is con- 
sistent what one would expect from observations of the IR lumi- 
nosity density and cosm ic star formation history, which have been 
found to peak at z ~ 1 (Hop kins & Beacorn]|2006l and references 
therein); at those redshifts the IR luminosity density is produced 
by roughly equal contributions from 'norma l' systems and LIRGs 
dLe Floc'h etal]|2005l ; IMagnelli et ai1l2009l) . The model further- 
more predicts that ULIRGs will become more important than 'nor- 
mal' systems in producing the CDIRB only at z > 2, and will 
equal that of LIRGs at even higher redshift z <; 3. However, the 
vast majority of the total CDIRB emission comes from objects a 
low to moderate redshifts, with ~ 50% arising from contributions 
at z <J 1 and 90% from z <J 2. Again, the discontinuities in these 
quantities at z ~ 2 are artifacts of the observed galaxy mass func- 
tions adopted, rather than a physical effect. 

In Figure [6] we isolate the contribution from sources above a 
given redshift threshold. This can then be compared to observa- 
tional estimates measured via a stacking analysis using BLAST 
data at 250, 350, and 500 pm which found roughly half of the 
CDIRB at those wavelengths was contributed by sources at z ^ 1 



6.2 Steady-state versus interaction-driven star formation 

The primary emission mechanism for producing the CDIRB is star 
formation, in which dust associated with periods and regions of ac- 
tive star formation reprocesses their UV-optical light into thermal 
emission. The model predictions presented in this work allow us to 
unfold the relative importance of two different modes of star for- 
mation: steady-state and interaction-driven (see § l2,ll and l2.2l for an 
overview). In Figure [8] we show the differential fraction of the to- 
tal CDIRB emission contributed as a function of redshift, including 
both the total and that arising from merger-driven bursts. We find 
that mergers generate only a small fraction, on the order of a few 
percent, of the emission at all redshifts. Steady-state star formation 
therefore contributes the vast majority of the CDIRB flux density 
observed today. The same is true at all wavelengths (see Figure|9]l, 
with merger-driven star formation accounting from only ~ 2 — 5% 
of the total emission. 

This owes trivially to the fact that, in our model merger- 
driven bursts account for a similarly small fraction of the IR 
luminosity density, which is consistent with a number of ob- 
servational constraints (out to z ~ 1; [ Brinchmann et al. 19981: 
i Bell et al.l2005l;|Lotz et alj2008bl;|jogee et al.ll2009l ; lRobaina etall 
2009; Lopez-Saniuan et al. 2010) and semi-empirical estimates 
(Hopk ins & Hernquistll2010l) . It is important to note, however, that 
in this context 'merger-driven' refers only to the star formation ow- 
ing directly to the interaction beyond what the two disks would 
have formed in isolation. When all on the star formation, includ- 
ing the steady-state mode, is taken into account, systems that 
would be identified as o ngoing mergers (e.g., via their morphology; 
lLotz et al .12008 a, 2009) can be a factor of 1.5-3 higher (also shown 
in Figures [8] and |9). However, even in this limiting cases mergers 
will never contribute more than ~ 10% of the CDIRB emission. 



6.3 The role of obscured AGN 

Though star formation clearly dominates the CDIRB, observa- 
tional constraints on the importance of obscured AGN remain a 
topic of some controversy and debate. For example, a number of 
authors dAlmaini et alj[l999l : iRisaliti et alj|2002l; iBlain & Phillips! 
l2002l ; lBargeretal.ll200ll) have claimed that AGN contribute of or- 
der 5-10% of the emission, making them comparable in importance 
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Figure 8. The differential contribution from merger-induced bursts (red) 
relative to the total (black) CDIRB as a function of redshift. It is impor- 
tant to note that in this case merger-driven refers to star formation directly 
attributable to the merger event, rather than all star formation in ongoing 
merger. The total contribution from ongoing mergers, including the steady- 
state mode is also shown (blue). 



Figure 10. The differential contribution to the CDIRB from obscured AGN 
relative to the total (black) for two different cases: a limiting case in which 
the total bolometric luminosity of all AGN is assumed to be reprocessed 
into the IR (blue), and the prediction from the fiducial model assuming the 
templat e SEP library and m odel parameters that best-describe Type II AGN 
Siebenmorgen et al. (2004b). 
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Figure 9. Same as Figure [8] but for the fractional contribution function of 
the observed CDIRB flux density as a function of wavelength. 



to merger-driven starbursts while others jSevergninie"tai]|2000l : 
ISilva et alJl2O0l ; lFardal et al]|2007l) claim they make a trivial con- 
tribution. Again, our model predictions can be used to separate out 
obscured AGN, and estimate their relative importance in producing 
the CDIRB. 

In Figure [10] we show the differential fraction of the total 
CDIRB as a function of redshift for the fiducial model (red shaded 
region) as well as a limiting case in which we assume that 100% 
of the AGN's bolometric luminosity is reprocessed into the IR at 
lO^im <, X <, 100/im. Therefore, our fiducial model predicts that 
AGN contribute 0.7 ± 0.2% of the total CDIRB emission, with a 
firm maximal contribution of 3 ± 1% in the limiting case. This 
results from a simple energetic argument: AGN do not produce 
enough photons to substantially affect the CDIRB. 

Though the model predicts the AGN contributes <> 1% of the 
total CDIRB emission, it is in principal possible that it could con- 
tribute relatively more at certain wavelengths. Therefore, in Fig- 
ure QT] we present the fraction contribution from obscured AGN to 
the CDIRB as a function of wavelength. In this case, rather than ex- 
amining the bolometric limiting case as we did previously (which 
by definition ignores wavelength dependence) we can instead vary- 
ing the model parameters about our fiducial choice (Av = 16 and 



R = 1 kpc) to examine their effect on the prediction. In the left 
panel, we hold R constant while varying Av, and in the right panel 
we hold Av constant while varying R. In most cases, the AGN con- 
tribution will peak at ~ 30 — 40^im owing to its relatively warmer 
IR SED as compared to star-forming systems of the same IR lumi- 
nosity (see examples in FigureQJ. Only in the most extreme column 
densities and largest sizes will the template SED be cold enough to 
contribute more at longer wavelengths. However, while the AGN 
can contribute somewhat more than ~ 2% at the shortest wave- 
lengths, they are clearly subdominant in general, even under the 
most extreme assumptions of obscuration and size. 



7 THE 7-RAY OPACITY OF THE UNIVERSE 

The CDIRB and diffuse backgrounds more generally also repre- 
sent the primary s ource of opacity to very high energy 7-rays in 
intergalactic space (Jelle vlll966l :l Fazio et al. 1966; Fazio & Stecker 
ll97fj[|Gould & Schrederlll967l : ISteckedll969l : IStecker et alj|l992l) 
When the mass-energy of photons is comparable to the rest-mass 
of electrons, they can interact via electron-positron pair production: 
7 + 7 — ¥ e + + e~ . The cross section for interactions between a 
photon with energy E-, with a background photon at energy e = hv 
at redshift z is: 
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where <tt is the Thompson cross section, e t h is the threshold en- 
ergy, n = cos 8 where 6 is the angle of incidence, and n e = 
dn-f/de is the differential comoving number density of background 
photons at z given by: 



e 2 n e (e, 2) = — vl v (v,z). 
c 



(12) 
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Figure 11. The fraction of the CDIRB flux density as a function of wavelength owing to obscured AGN for the full range of SED model parameters. The left 
panel shows the contribution at fixed R = 1 kpc for (bottom to top) Ay = 1, 2, 4, . . . 128 ( log(AT H /cm~ 2 ) = 22.3, 22.6, 22.9, . . . 23.5). The right 
panel shows the contribution at fixed A v = 16 (\og(N H /cm -2 ) = 22.5) for R = 0.25, 0.5, 1.0, 16 kpc. 



To calculate the optical depth of a photon with energy in 
the observers frame that was emitted by a source at redshift z em , 
simply integrates along the line of sight: 



. , dt 
az c— 

dz 



-l 1 — //' f°° 

dfj,' — - — / de'n e (e' , z)a~ ll (E 1 ,e , jx) 



(13) 



where cdt/dz is the comoving line-of-sight distance (see Hogg 
1999). Since the primary emission mechanism at these energies 
is thought to be synchrotron self-Compton, and therefore is 
well-d escribed by a simple power-law (e.g., lAharonian et al.l 
l2002bl) . one can used <r 77 inferred from the observed GeV-TeV 
SEDs of 7-ray s ources to constrain the intensity of the dif- 
fuse background dStecker et alj 1 19921 : Istecker & de JageJ fl993l; 
Steckeret alJ l200l iDwek & Slavinl Il994l: |Pwek & Krennricbl 
20051 : Ide Jager et all 119941: iMadau & Phinneyl Il996l: iFunk et all 
1998l:|Primack et alJl999l:lBullock et al.ll999l:lRenault et alj200ll; 
Aharonian et al. 2002b, 2006; Schroedter 2005; Katarzynski et al. 
20061: iMazin & Rauej 120071; iGeorganopoulos et alj 120081 ; 



Finke & Razzaque 2009). In general, for a uniform background of 



photons, the cross section peaks at /3 = 0.7 which corresponds 
to a background photon with wavelength E 1 e ~ 4(m e c 2 ) 2 or 
(A e //zm) ?s 1.2 x (S 7 /TeV). The diffuse background has a 
trough in the near-IR b e fore rising ste eply out to A 6s ~ 200/^m 
dHauser & Dwekl 1200 ll ; lLagache et alj 120051) . Therefore, the the 
opacity of the Universe to ~few x TeV photons (typically referred 
to Very High Energy, or VHE 7-rays) is generically dominated by 
the mid-IR component of the diffuse background. 

Since our model predicts the CDIRB intensity, and therefore 
number density, we can invert this procedure and instead provide 
predictions for opacity of the Universe due t o pair creation inter- 
actions as a function of E^ and z em (see also Primac k~et alj|2005l 



2008; Franceschim et alj2008l : lFinke etail2 009). In Figure[l2] we 
compare the model predictio ns to estimates of t 77 usin g the 'av- 
erage' local CDIRB model of IDwek & Krennricbl J2005h . We find 
good agreement between the two, and also confirm that by only in- 
cluding the IR component of the diffuse background emission, we 
are essentially capturing the opacity of the Universe to observed 
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Figure 12. The opacity of the local Universe to electron-positron pair pro- 
duction interactions (7 + 7 — > e + + e~) as a function of the observed 
energy (Ej), assu ming a source at z = 0.0 3. The results assuming the 
mean CDIRB from Dwek & Kr ennrichl fc005l) are shown as a solid black 
line, and the same only including the optical through near-IR component 
of the background photons (A < lOfim) as a dashed line. We also show 
the model predictions (red shaded region), which includes background pho- 
tons with A > lOfim. Over the range where the CDIRB dominates the 
7 + 7 opacity (Ej > 1 TeV), the model provides a good match to the 
IDwek & Kre nnrich (2005) results. 



photons with _E 7 J> 2 — 3 TeV; below these energies the model 
predictions should be considered a lower limit. 

Figure [13] shows the absorption edges predicted by the model 
for a number of optical depths ranging from r 77 = 1/3 — 100. Two 
TeV sources have b een observed to have exponent ial cut-offs in the 
VHE 7-ray SEDs dAharonian et alJll999L f2002ah which are con- 
sistent with the model predictions for r 77 <; 1 at those redshifts. 
These results suggest that the Universe is very optically thick to 
TeV photons beyond a redshift of z em ~ 0.2 with r 77 > 3 for 
£ 7 ^ 3 TeV and r 77 Z 10 for E 1 Z 8 TeV. Furthermore, at 
Zem 0.5 the Universe is optically thick down to ~ 1 TeV, and 
r 77 w 10 down to ~ 3 TeV. At higher energies (_E 7 <; 10 TeV) 
the Universe is very optically thick (r 77 ~ 3 all the way down to 
z ern ~ 0.01. Though significant improvements in sensitivity, par- 
ticularly at <; 10 TeV, are planned in t he next generation of ICATs 
(e.g., the Cherenkov Telescope Array; iHermann et al ] |2008h . these 
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Figure 13. Model predictions for the absorption edges due to electron- 
positron pair production as a function of redshift or line-of-sight comoving 
distance (left/right absciss) and observed photon energy (E~ ) for a range 
of optical depths: r 77 = i, 1, 3, 10, 30, and 100. We also include the 
cut-off energies observed for two local TeV sources jAharonian et al.l 19991 
l2002ah . 



results suggest that the high optical depth to pair creation might 
pose a significant problem for VHE observations at even somewhat 
low redshifts. One can, however, expect to observed ~ TeV photons 
out to fa 100 - 150 Mpc. 



8 CONCLUSIONS 

We present a physical model for the origin of the C DIRB utiliz- 
ing the semi-empirical framework of Hop kins et all J2010bl) . The 
model tracks three distinct sources of IR emission: steady-state star 
formation, interaction-induced starbursts, and obscured AGN activ- 
ity. We also include all the relevant systematic uncertainties, which 
are dominated by the dynamic range in observational estimates 
of the stellar mass function (and increasing with redshift). The 
IRLFs generated by this model, combined with a library of tem- 
plate SEDs for starb ursts dDale et alJkOQll ; [Chary & Elbazl[200l| ; 



Lagache et al 



1993 



2001 



20031) and obscured AGN ( Siebenmorgen 1991, 
Siebenmorgen &~Kruegel| [l992; Siebenmo rgen et al]|l992l 
l2004allbl) , provide an excellent match to observations of the 



CDIRB from X b s ~ 10 — 1000/im. In contrast of alternative tech- 
nique s - inc l uding backwards ev o lution (e.g., iRowan-Robinsonl 



200 ll 120091: lLagache et all l200l |2004 120051; |Xu et all 120031 ; 
Franceschini et alj 2008 ; Finke & Razzaque 2009l). semi-analyti c 



l Primack et al. II 19991 120051. 120081: iDevriendt & GuiderdonfeOOOl) . 
and CCE dPei & Fallll 19951 : |Pei et al.lll999l) models - our approach 
provides a robust context in which to examine the physics its pro- 
duction. 

By varying the model assumptions, we find that the normaliza- 
tion and peak of the CDIRB emission are determined in large part 
by the strong redshift evolution of the gas content of steady-state 
star-forming galaxies. Without this evolution, the model vastly un- 
derpredicts the total intensity of CDIRB and peaks at shorter wave- 
lengths than is observed, owing to the lack of a significant con- 
tribution form high-luminosity and high-redshift disks. Therefore, 
the observed CDIRB can be taken to reflect this strong evolution 
in the gas content of disks - an effect whic h has been seen in nu- 
merous samples of individual objects (e.g.jBell & de Jondl200ll; 
Kannappaj2004lMcGaughl2005klshaplev et al.l2005l ; lDaddi et al.1 
20091 : iTacconi et alj2010l) 



We can also use the model predictions to determine the relative 
importance of IR-luminous galaxies in producing the CDIRB as 
a function of redshift, luminosity, and emission mechanism. The 
model results indicate the following: 

(i) The CDIRB is primarily produced by equal contributions 
form objects at z ~ 0.5 — 1 and z 1, in agreement with re- 
cent observations by the BLAST experiment. However, in the ob- 
server's frame wavelengths past the peak contain an larger contri- 
bution from z <; 1 than those at shorter wavelengths. 

(ii) Most of the CDIRB is contributed by normal galaxies and 
LIRGs, though the ULIRG contribution becomes significant at 
high-redshift z Jj 2. This also leads to a larger contribution from 
ULIRGs at wavelengths past the peak. 

(iii) Ongoing mergers contribute less than 10% of the total emis- 
sion at all wavelengths, with less than ~ 1 — 3% produced by the 
merger-driven burst itself. 

(iv) Obscured AGN account for < 1% of the CDIRB at all 
wavelengths, with a strong upper limit at 4% of the total emission. 

In the future, these predictions can be tested via deep, wide sur- 
veys with the Herschel Space Telescope (e.g., the Herschel Multi- 
tiered Extragalactic Survey, or HerMESf] combined with stack- 
ing analyses similar t o those presented bv lDevlin et alj d2009l) and 
iMarsden etakl d2009l) . 

The CDIRB also represents the primary source of opacity 
for VHE 7-rays, owing to electr on-positron pair production in 



teractions (7 + 7 — >■ e + + e~: I Jellevl 1 1 9661: iFazio et al l 
Fazio & Steckeil Il970l : iGould & SchrederJ 1 19671: ISteckerJ 



1966 



1969 



Stecker et al.lll992TT Thus, our results also provide predictions for 



the opacity of the Universe to the most energetic photons. We find 
that the model predictions ar e consistent with high-energy cutoffs 
for TeV sources at z ~ 0.03 dAharonian et alJI 19991 [2002a). They 
also indicate that while the Universe is highly opaque to E-, <; 10 
TeV photons (t 77 J> 10 at z ~ 0.06), with the next generation of 
ICATs we can reasonably expect to detect sources out to ~ 50—150 
Mpc. 
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